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Abstract. We introduce a class of dynamical systems having an invariant 
measure, the modifications of well known systems on Lie groups: LR and 
L+R systems. As an example, we study modified Veselova nonholonomic rigid 
body problem, considered as a dynamical system on the product of the Lie 
algebra so{n) with the Stiefel variety Vn,r, as well as the associated eL+R 
system on so(n) X Vn,r- In the 3—dimensional case, these systems model the 
nonholonomic problems of a motion of a ball and a rubber ball over a fixed 
sphere. 


1. Introduction 


This paper describes a class of dynamical systems allowing an invariant measure, 
a modification of LR systems introduced by Veselov and Veselova [2T] and L+R 
systems introduced by Fedorov (see mm)- In particular, they model the non¬ 
holonomic problems of rolling the Chaplygin ball and the rubber Chaplygin ball 
over a spherical surface. 

Recall that the motion of the Chaplygin (balanced, dynamically asymmetric) 
ball over a fixed spherical surface is described by the equations 


( 1 ) 


dt 


k X : 


|7 = e7X+, 


where oj is the angular velocity of the ball, 7 is the unit vector directed from the 
center of the fixed sphere to the point of contact, k = Iw + Doj — D{uj,^)j is the 
momentum of the ball with respect the contact point, and I = diag(/i, I2,13) is the 
inertia operator of the ball (e.g., see laiz])- 

The parameters D and e are equal to rri^p and aj{azLp), where m is the mass 
of the rolling ball, p its radius and a is the radius of the fixed sphere. The sign 
denotes the rolling over outer surface of the fixed sphere, while the sign ” 
denotes either the rolling over inner surface of the fixed sphere {a > p), or the case 
where the fixed sphere is within the rolling ball (a < p, in this case the rolling ball 
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is actually a spherical shell). As a tends to infinity, e tends to 1 , and we obtain the 
equation of the rolling of the Chaplygin ball over a horizontal plane. 

In the space (w, 7 ), the density /r of an invariant measure is equal to 

(2) = v'det(I + D¥.) (1 - Z?( 7 , (I + 


the expression given by Chaplygin for e = 1 [5], and by Yaroshchuk for e 7 ^ 1 [^ . 
Through the paper the operator E denotes the identity operator on the appropriate 
spaces. 

Similarly, the motion of the rubber Chaplygin’s ball over a fixed spherical surface 
is described by the equations (Ehlers and Koiller [11] , see also Borisov and Mamaev 

00 ) 


(3) 


d d-. „ 

—m = mxa; + A 7 , —7 = 67 x 0 ;, (o;, 7 ) = 0. 


where m = (I+£>E)o; = Iw, A = (m, I~^ 7 )/( 7 , 1 “^ 7 ). Here the constraint (w, 7 ) = 0 
model the ’’rubber” property of the ball: the rotations of the ball around the normal 
to the spherical surface at the contact point are forbidden. 

Note that the above system is well defined on the whole phase space (to, 7 ), and 
(a;, 7 ) = const denotes its first integral. In the space (to, 7 ), or (a;, 7 ), the system 
has an invariant measure with the density given by 


(4) 


Me = (I ^7,7)’^C 


see [mno] for 6 = 1, and m for 6 7 ^ 1. 

The existence of an invariant measure for some nonholonomic systems is an 
important property related to the geometry of the problem, its Hamiltonization, as 
well as to the possibility to solve the problem by quadratures (e.g., see [niiiiini 

unidiiioiiia]). 

We shall prove that the densities ([2]) and (|4|) are particular cases of the densities of 
invariant measures of 6 -modified L+R and LR systems, respectively (see Theorems 
1, 3 and 5, and Examples 1 and 2). Eurther, as a specific example, we give the 
expression of an invariant measure for the modified Veselova nonholonomic rigid 
body problem (Theorem 4), considered as a dynamical system on the product of 
the Lie algebra so(n) with the Stiefel variety 14i,r, as well as the expression of an 
invariant measure for the associated cL+R system on so(n) x Vn^r (Theorem 6 ). 
Eor r = 1, these systems represent natural multidimensional generalizations of the 
equations m and 


2. LR AND L+R SYSTEMS 

2.1. LR systems. Let G be a n-dimensional compact connected Lie group G, 
g = Lie{G) its Lie algebra, (•,•) an Adc-invariant scalar product on g, and let 
(•, -(i be a left-invariant metric on G given by the positive definite operator (called 
inertia operator) I : g g ^ g*: (?7i,??2)i = (I( 5 “Ni),e TgG. Here 
we identify g and g* by (•,•). 
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The LR system on G is a nonholonomic Lagrangian system (G, L,X>), where 
L = is a left-invariant Lagrangian and I? is a right-invariant nonintegrable 

distribution on the tangent bundle TG, determined by its restriction 0 to the Lie 
algebra [21]. Let 1) be the orthogonal complement of h with respect to (•,•). Then 
the right-invariant constraints, in the left-trivialization of TG, can be written as 
w £ Adg-i 0, or (w, Adg-i f)) = 0, where w = g~^ ■ g is the angular velocity. 

The LR system (G, T, T) is described by the Euler-Poincare-Chetayev equations 
on T*G{m,g) (or TG{u},g)) 

k 

(5) m = [m,w]-b ^ A'* Cs , g = guj, 

S^l 

where m = Iw S g* = g is the angular momentum, = g~^ - Eg ■ g, Ei,... ,Ek 
is the orthonormal base of 1), and A® are Lagrange multipliers which can be found 
by differentiating the constraints 4>s = (w, Cg) = 0, s = 1,..., k. These equations 
define a dynamical system on the whole cotangent bundle T*G, and the right- 
invariant constraint functions (ps are its first integrals. Also, to the system ([5|) we 
can associate a closed system on g™+^(TO, ei,..., Cfe): 

k 

(6) m = [m,a;] 4-^A^es, Cs = [cs^uj]. 

s^l 

The LR system (jH]) possesses an invariant measure g dm A dei A • • • A de^ with 
density (see m) 

(7) /r = •\/det(eg,I“ie;), s, 1 = 1,...,/c, 

implying that the original system ([5]) has an invariant measure y^det(I“^|g-i[,g) 17". 
Here is the restriction of the inverse inertia tensor to the linear space 

g~^tig C g, by dm, de^ we denoted the standard measures on g with respect to the 
metric (•, •), and 17 is the standard symplectic form on T*G. 

2.2. L+R systems. In addition to the inertia operator I defining the left-invariant 
metric (•,-)ii introduce a constant linear operator n° ; g —g defining a right- 
invariant scalar product (•, -(n on G: 

(?7i,??2)n = {Lfr]ig~^,r]2g~^) = {Ugf]ig~\g2g~^), Hg = Adg-i on° o Adg, 

Vi^V 2 G TgG. We suppose that Kg = I + Hg is nondegenerate and positive definite 
on the whole group G. The L+R system on G is defined as a dynamical system 

(8) ^ (Iw + Hgo;) = [Iw + Hgo;,!*;], g = g-w. 

This is the the modification of the geodesic motion on the group G with respect to 
the metric (•, •)! + (•, •)n, by rejecting the term [a;,ngw] at the right hand side of 
the first equation in (|S|). 

In the view of the definition of Hg, its evolution is given by the matrix equation 
Hg = Hg o ado; + ad^ oRg. Note that for compact group we have ad^ = — ad;^, and 
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therefore, we get a closed system 

(9) ^ (Iw + nw) = [Iw + IIw, w], 11= [n, adij] 

on g X Sym(g), where Sym(g) is the space of of symmetric operators on g, which we 
also refer as a L+R system. It possesses the kinetic energy integral i(Ia; + Ho;, w) 
and an invariant measure /r dw A dll with density 

^ = \/det(I + n), 

where dll is the standard measure on Sym(g) [14]. It appears that every L+R 
system can be seen as a reduced system of a certain LR system on a direct product 
G X g", where g is considered as a commutative group (see Theorem 3.3 in [18jL 

3. Modified LR systems 

As the sphere-sphere problems dTj), dSj) suggest, it is natural to consider modifi¬ 
cations of the equations m and We define the eLR system on the space 

g''+^(a;,ei,... ,efc), or g'=+^(TO, Ci,..., Cfc), 

by the equations 

k 

(10) th = [to,w] + A, A = ^A*ei, m = Iw, 

i=l 

(11) ii = e[ei,a;], 

i = 1,... ,k. The term A can be interpreted as a reaction force and it is determined 
from the condition that the trajectories (w(t), ei(t),..., ek(t)) satisfy constraints 

(12) (j)i = (w, ei) = Ci = const, i = 1,... ,k. 

By differentiating (TT^ . we obtain the linear system 

(13) e([ei,a;],a;) + (ei,I“^[m,a;]) = 0, A^-= (ci,I'^e^), 

3 

and we get the Lagrange multipliers in the form 

(14) A* = -^(e„I-i[m,a;])A^A 

3 

Here A*1 is the inverse of the matrix Kij = (I“^ei, Cj), i,j = l,...,k. In particular, 
A does not depend on e. 

Theorem 1. The eLR system (nni), (HD), (HI has an invariant measure 

(15) /ie dm A dei A • • • A de/c, i.e., /i^ dw A dei A • • • A de^ 


where the density is given by 

(16) ^e = (A)i, A = det(Aij) = det((I“^ej,ej)), 


i,j = l,...,k. 
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Proof. Consider the vector field 

X = (to, ei,... ,efe) = {[m,uj] + A, e[ei, w],..., e[efe, a;]). 

The volume form dTSl) is invariant with respect to the flow of (flOl) if and only if 

Cx dm A dei A • • • A de^ = 0, 


i.e., if fie satisfies the Liouville equation 


(17) fie (div(TO) + div(ei) H-h div(efe)) + /i^ = 0, 


where div(-) is the standard divergence on g. It is clear that div(ei) = 0,i = 1,... ,k. 
In Theorem 1 [21], it is proved that 

k 

(18) div(TO) = div (A) = — ^ A*l([I“^ei, I“^to], e^). 

i,i=i 


For any regular matrix A we have the identity 


(19) 

Thus, 


— det A = det Atr(A ^A). 



ij'=l 


= ^ ((I [ej,uj]) + (I [ei.uj])) 

k 

= Me X! 


which together with (fT51) implies (I17p . Since dm = detl-dw = const-dui, if instead 
of TO we take the variable w, the expression for an invariant measure remains the 
same. □ 


3.1. Momentum equation. As above, let Ei,... ,Ek be an orthonormal base of 
(). Further, let E^+i ,..., En be an orthonormal base of i) and OEt be the adjoint 
orbit oi Ei, i = 1,... ,n. It is clear that we can also consider the eLR system (US, 
(HD), (HD) on the space 

M = {{uj,ei ,... ,e„) G g X Oe^ x ■ x Oe^ \ (e*,ej) = %}, 
simply by taking i = I,..., n in the equation (HU. Then it has an invariant measure 
(20) /Te dw A dei A • • • A dcnlvvi- 

Moreover, on Ai we have well defined orthogonal projections pr^, pr^j (pr^ + pr^j = 
E) onto the complementary subspaces of g: 

T-L = spanjei,..., Cfc}, V = span{e/c+i,..., e„}. 
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which, according to (HU, satisfy the equations 

(21) = e[pr 2 ,,adc^], ^ pr„ = e[pr„, ad^]. 

Note that with the above notation, we can express A in (flUl) as 

A = — pr-^ w], 

where : TL ^ TL is the inverse of the mapping A = 1“^ j-H = pr-j^ : TL TL. 
Following m, let us introduce the momentum 

m = pr^, luj + pr^ lu = Jlu, J = pr^, I + pr^ = E + pr 25 (I — E). 

Consider the space 

J\f {(iTli ^fc+l 5 ■ • ■ ) ^n) ^ 0 ^ ^Ek+i X • • • X Oj^ri I ) ^ij } ■ 

Proposition 2. The eLR system on M has the following form 

(22) m = e[m,a;] + (1 - e) pr 25 [Iw,u;], 

ei = e[ei,uj], i = fc + l,...,n. 

Proof. From the equations m and (EU, the evolution of the momentum m reads 

rh = ^ (pr-D + pr^i w) 

= e (pr 2 ,[a;,Icj] — [w, pr^ Icc]) + pr^, ([Iw, w] + A) 

+e(pr„[a;,a;] - [a;,pr„a;]) +pr„w 
= e[m, w] + (1 — e) prjjpo;, w] + pr^ ch . 

On the other hand, from m it follows: e{[ei,a;],a;) + {ei,uj} = {ei,Lu) = 0, 
i = 1,..., /c. Whence 

k 

(23) pr^ uj = ^(w, ei)ei = 0. 

i=l 

Finally note, since the linear operator J is invertible, ui = J~^m and we have a 
closed system (l22l) on Af. □ 

In the case e = 1, the first equation in (1221) reduces to the momentum equation 
(2.7) of [15]. Also, as in Theorem 3.2 [T5], we have 

(24) det I • det A = det I • det(I“^ |^^) = det J = det(I|i 5 ), 
where Iji? = pr^, ol : V ^ V. Therefore: 

/Te dw = /ie det 7 -— dm 
om 

= (det(I“^|?i)) (det J)“^dm 

= (det 1)“^ (det I|x))^“^dm . 

Combining the expression of the invariant measure (1^01) on the space M and the 
above identity, we get the following statement: 
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Theorem 3. The eLR system (12211 has an invariant measure 

/ie dm A dcfc+i A • • • A de„|A7, 
where the density is given by 

fie = (detl|r>)^"^ = (det({Iei,ej)))A“\ i,j = k + l,...,n. 

Example 1. Let us consider the case when TL is the isotropy algebra of 7 = ei: 
H = {x € g I [x, 7 ] =0}. Then V can be identified with the tangent plane T^O of 
the adjoint orbit O through 7 . Since TL and V are uniquely determined by 7 , we 
can write the closed system in variables (m, 7 ) or ( 01 , 7 ): 

(25) rh = e[m,w] + (1 - e)pri,[Ia;,a;], 7 = e[ 7 ,w], 

m = prj, Iw + pr^ ui. 

Also, in the special case when Tl is one dimensional, spanned by 7 = Ci, we have 
the equations 

(26) m = [m,uj] + X'j, 7 = e[ 7 ,a;], 

A =-(7,I-i[m,w])/(I“i7,7), to = Iw. 

The above examples coincide in the case of the Lie algebra so(3). Under the 
usual isomorphism between the Euclidian space and so(3) 

/ 0 -A3 A2 \ 

(27) A = (Ai,A 2 ,A 3 )^A= A 3 0 -Ai , 

\-A2 Ai 0 J 

replacing the inertia operator I by I = I + DE, the equations (|26l) recover the 
equations © of the rubber Chaplygin ball, while the expression for the density of 
the invariant measure (HU) recovers the density Also, according to (1^ . we can 
write the rubber Chaplygin ball equations in the equivalent form 

(28) = ern x w + (1 — e) (Iw x w — (Iw x w, 7 ) 7 ), -^7 = 67 x w, 

dt dt 

rn = Itj + ( 7 , a; — Iw) 7 . 

3.2. Modified Veselova problem. Following [15], we define e-modified Veselova 
nonholonomic rigid body problem as follows. Let ei,...,e„ be a (moving) or¬ 
thonormal frame of the Euclidean space (M", (•,•))• Consider the orthogonal de¬ 
composition 

(29) soln) ='Hr® TTr, 

Tir = spanjep /\eq\r < p < q < n}, Vr = spanje^ A e_,- 11 < i < r, 1 < j < n}. 
Then 

WVrV = rf&soin), 

where 

r = ei (g) ei -I- • ■ • -b e^ Cl Bj.. 
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The above projection depends only on the point (ei,..., e,.) of the Stiefel variety 
Vn^r, realized as a submanifold of ]R"'’'(ei,..., e^) by the constraints 

Vn,r ■ 1 ^ — bj — 

Therefore, instead of using the variables m,ei A e^ , 1 < i < r, 1 < j < n, we can 
write the equations (l 22 l) on 


so{n) X Vn,r (m, ei, ...,er) 


as follows: 

til = e[m,a;] + (1 — e) (r[Iw,a;] + [Iw, wjT — r[Iw,a;]r), 
(30) Bi =-euJBi, i = l,...,r, 

m = w + r(Iw — w) + (Iw — a;)r — r(Iw — a;)r. 


Along the trajectory (m(t), eier(t)) of the modified Veselova problem 
(Pl) . the angular velocity matrix uj{t) has the form 


uj(t) = 


I 0 


-iOxrii) 




OJirit) 




^rn{t') 


Cr 


^rn {i') 




where uJij(t) = {ei{t) A ej{t),uj(t)) and Cr is a constant {n — r) x (n — r) matrix. 
Next, define the left-invariant metric via the relation 


(31) I(Ei A Ej) = UiUjEi A Ej, 

where Ei,..., E„ is the standard orthonormal base of M" and A = diag(ai,..., a„) 
is positive definite. 

Then, from Theorem [3] and Theorem 5.1 of m, we get: 

Theorem 4. The e-modified Veselova system dSO]), dSlI) has an invariant measure 

■ • - 0^(61 A ■ • • A Brfi J dm A dei A ■ • ■ A derU(„)xy„,,, 

where the summation is over all r-tuples I = {1 < ii <■■■< C < n}, and 
(ei A • • • A er)i are the Pliicker coordinates of the r-form ei A • • • A e^. In the case 
r = 1, the density is proportional to (ei, Aei)*^A“i)("“ 2 )_ 


Here dei A • • • A de^ is the standard volume form on K”’'. Note that, for the 
inertia operator I replaced by I = I -I- HE, r = 1 and n = 3, the equations (1501) give 
another form of the rubber Chaplygin ball equations (l28ll . 
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4. Modified L+R systems 
The eL+R system on the space 

0 X Sym(0)(a;,n) or 0 x Sym(0)(k, R) 

is defined by 

(32) k=[k, w], n = e[n, ad(^], k = Iw + Ho;. 
Theorem 5 . The eL+R system possesses an invariant measure 

/rdwAdn i.e., dk A dll, 

with density same as in the case of the usual L+R systems: 

(33) fi = •\/det(I + n), 
where dll is the standard measure on Sym(0). 


Proof. The proof is a modification of the corresponding statement for L+R systems 
(see [Mills]). We have 

k = (I + n)d; + nw = (I + n)a; + en[a;, w] — e[uj^ Rw], 

Thus, we can represent the equations (1321) in the form 

(34) w = (I + R)“^ ([Iw,a;] + (1 - e)[Ra;,a;]), R = e[R, ad^;]. 


The volume form ((551) is invariant with respect to the flow of (1551) if and only if 

(35) pL ^div(+) + div(R)^ + /i = 0, 

where we take the standard divergence on 0 and Sym(0): 

i 

Rere we take coordinates of w and R with respect to the orthonormal base Ei,, En 
of 0 and the associated base Ei (g) Ej of linear operators on 0. 

It is clear that div(R) = 0. Define n x n matrix R: 

^ _ 5((I + R)w)i _ 5([Ia;,a;] + (1 - e)[Rw, w])i , i 

- r) - r) 'I - i, . . . , 71. 

OUJj OUJj 

Then 

R = - ad(^ oI + aditj +(1 - e)(- ad,^ oR + adn^j) 


and we can write 

div(d;) = tr((I + R)“^R). 

In view of symmetry of I + R, the skew symmetric part of R does not contribute 
to the expression for the divergence. 

Taking into account dSH), the symmetric part of A has the form 

R+ = ^ (I ° s-do; — ad,^ ol + (1 — e)(R o ad,^ — ad^^ oR)) 

= i ^(I + R) o ad^ - ad^ o(I + R) - R^ . 
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As a result, we obtain 

div(a;) = tr((I + n)“^n-|.) 

= i tr ^(I + n)“^((I + n) o ad(^ - ad(^ o(I + 11) - II)^ 

= -itr((i + n)-iri) =-itr(^(i + n)-i|(i + n)) 

= —-(det(I + n))~^ — det(I + n) 

= -(x/det)! + n))-!^ Vdet(I + H), 
at 

where we used the unimodularity condition for compact groups tr ad^^ = 0 and the 
well-known identity (HU. We conclude that fi = -\/det(I + II) satisfies the Liouville 
equation (l35)l which establishes the theorem. □ 


Remark 1. Note that the kinetic energy iJ = i(k, oj) is conserved for eL-|-R systems 

— (k, WjO;) = 2 ((I -I-n)w,a;)-I-( — (I -b n)a;,a;) 

= 2([Ia;, to] -b (1 — e)[nto, to], w) -b e(n[to, to] — [to, Ilto], to) = 0, 

while, for the eLR systems, the kinetic energy H = ^(Ia;,a;) is conserved only on 
the invariant submanifold (j)s = (w,es) = 0, i.e., when pr-^co = 0. However in the 
case e = 1 , we have also the preservation of the following modification of the kinetic 
energy: 

F = - (pr.„to,Ico). 

Indeed, -^H = {to. A), and from (1^ . (l23l) . we have 

^(pr„a;,Ia;) = (epr.„[to, to] - e[a;, pr.„ to], Ico)-b (pr^ to, [Ito, to]-b A) 

= (1 — e)(pr.^ CO, [Ito, to])-b (to. A) = (to, A) for e = 1. 


Remark 2. The (e-modified) LR and L-bR systems can be considered on non¬ 
compact groups, when we have an invariant measure for unimodular groups as 
well. Recall that the group G is unimodular if tr ad^j = 0. 


Example 2. As in Example 1, consider the orthogonal decomposition g = ©X>, 

where H is the isotropy algebra of 7 = ei. Let us take H = Dpr-p. Then the 
modified L-bR system (l32|) takes the form 

k=[k, to], 7 = e[7,to], k = Ico-b I^pr^, to. 

After the identification so(3) = given by (P71) . Theorem 3 recovers the invariant 
measure ©■ Another natural choice for the operator H is 

n = L>[[7 ,w],7] 

(see [mus]), which also leads, for g = so(3), to the sphere-sphere problem ([T]). 
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4.1. eL+R system on so{n) x 14,,r- Here we use the notation of Subsection 13.21 
Consider the decomposition of so{n) given by (1^ . and take H = Upr^^. As a 
result, we obtain the modified L+R system on the space so(n) x V4,r 

k = [k,w], 

(36) dj = -ewej, i = 1,... ,r, 

k = Icti + -D(T'oj + cuL — Lcur), L = e^^ (S) + ■ • • + e^ 0 

For r = 1 and e = 1, the equations (1361) model the problem of rolling without 
slipping of a Chaplygin ball over the hyperplane in R", orthogonal to ei (see [T^ 
[H]). As in |19j (see eq. (49) of [19)1. consider the metric defined by inertia operator 

(37) I(E,AE,-)= E,AE„ 

D - aittj 

where 0 < GiUj < D, i, j = 1,... ,n. 


Theorem 6. The eL+R system (1361) . with the inertia operator given by dSZl) has 
an invariant measure 



(ei A • ■ • A er)f 

* ’ ' ^ir 


-i(n-r-l) 

dk A dei A • ■ • A der|so(«)xV„.r- 


For r = 1, the density is proportional to (ei, A ^ei) 2 ^" 


Proof. Motivated by the relationship between 3-dimensional Veselova problem and 
the Chaplygin ball problem established by Fedorov m, define the operator I and 
matrixes w, m by 

I = E + I1I“^, w = Iw, m = prj, Iw + pr,j^ w. 


Then 

m = pr^, loj + prj, + pr.^ lui = Iuj + D pr^ w = k. 

Therefore, by using (1^ . we get 

(9k 9ni d'w 

det(I + n) = det — = det —— • det = det(I|-p^) • det I. 

OU! OW OU} 

Now, the definition of I can be seen as follows: it implies the identity I(E,iAEj) = 
Da~^aJ^EiAEj. By combining the above expressions with Theorem[S]and Theorem 
5.1 of [15] we obtain the statement. □ 


Remark 3. The system m is integrable by the Euler-Jacobi theorem [T] for e = 1 
(Chaplygin [9], see also [TJ |B]) and for e = —1 (Borisov and Fedorov |3], see also 
mm)- Similarly, in the case of the rubber rolling, we have integrability for e = 1 
(see mi uni) and e = — 1 (see 13 [S]). The problem with the addition of potential 
forces is studied in details in |8]. On the other hand, the integrable models of the 
rolling of the Chaplygin ball (on the zero level set of the SOin — 1)-momentum 
map) and the rubber Chaplygin ball over a hyperplane in R" are given in [19] and 
[13 ng, respectively. It would be interesting to study appropriate problems where 
we have a spherical surface instead of a hyperplane. The natural n-dimensional 
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variants of the equation o and ([3|) are, respectively, the equations ([36)l and dSoj 
(with I replaces by I = I + £)E), where we set r = 1. 
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search was supported by the Serbian Ministry of Education and Science Project 
174020 Geometry and Topology of Manifolds, Classical Mechanics, and Integrable 
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